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The rotation frequencies of young pulsars are systematically below their theoretical Kepler limit. 
R-modes have been suggested as a possible explanation for this observation. With the help of semi- 
analytic expressions that make it possible to assess the uncertainties of the r-mode scenario due to 
the impact of uncertainties in underlying microphysics, we perform a quantitative analysis of the 
the spin-down and the emitted gravitational waves of young pulsars. We find that the frequency to 
which r-modes spin down a young neutron star as well as the characteristic gravitational wave strain 
amplitude are extremely insensitive both to the microscopic details and the saturation amplitude. 
Comparing our result to astrophysical data, we show that for a range of saturation amplitudes r- 
modes provide a viable spindown scenario and that all observed young pulsars are very likely already 
outside the r-mode instability region. Taking into account the finite observation time, we find that 
the signal to noise ratio for gravitational waves is smaller than previous estimates, but large enough 
to detect such sources with next generation detectors. 



I. INTRODUCTION 

A neutron star is a complex system whose behavior is 
influenced by physics on scales from the Fermi scale to its 
size and whose description involves all forces of nature. 
This holds in particular when dynamical aspects involv- 
ing mechanical oscillations coupled to radiation fields are 
considered and a description should thereby intricately 
depend on all these details. Yet. the surprising find- 
ing presented in this work is that for some of the most 
relevant observables analytic expressions can be derived 
that are strikingly insensitive to these complicated mi- 
croscopic as well as macroscopic details. 

The rotation frequencies of pulsars, as well as their 
time derivatives, are among the most precise observa- 
tions in physics. These frequencies change over time 
due to magnetic braking and other possible spindown 
mechanisms like gravitational wave emission due to de- 
formations or oscillation modes of the star. In particular 
the unstable r-modes [TJ |5] are interesting because they 
must be stabilized by a viscous dissipation mechanism 
and thereby can probe the microphysics deep inside the 
star. It has already been shown |31 H] that important as- 
pects of the instability regions of r-modes are very insen- 
sitive to quantitative microscopic details but do depend 
on qualitative differences of the various possible forms of 
dense matter. Here we extend this study to the dynam- 
ical r-mode evolution of neutron stars and show that in 
this case the insensitivity to unknown input parameters 
is even more pronounced. In particular we derive semi- 
analytic expressions for the final frequency of the spin- 
down evolution and the spindown time of young neutron 
stars and show that these quantities are likewise insensi- 
tive to the particular r-mode saturation mechanism. 

Magnetic dipole radiation is considered as a standard 
mechanism for the spindown of pulsars and is used to 
determine an approximate characteristic spindown age 
[5]. In contrast, in this work we will instead only con- 
sider the spindown torque due to r-modes, which is only 
non-vanishing as long as the star is within the unstable 



region. Our results therefore provide upper limits on the 
actual frequencies of observed pulsars and should apply 
as long as the r-mode spindown dominates. However, 
they might not be applicable to magnetars where this is 
not guaranteed. A more detailed analysis which includes 
both electromagnetic and gravitational radiation [6J and 
thereby allows to make contact to pulsar ages will be 
given elsewhere. 

The fastest young pulsar observed so far is PSR J0537- 
6910 with a rotational frequency of /«62 Hz, which has 
been associated to the remnant N157B of a supernova 
that exploded in the Large Magellanic Cloud • We find 
that this spin frequency is just below the final frequency 
due to pure r-mode emission of a neutron star with stan- 
dard modified Urea beta equilibration processes, so that 
r-mode emission could indeed provide a quantitative ex- 
planation for the low spin frequencies of young neutron 
stars. If the r-mode scenario is realized, the fact that all 
observed stars are below the r-mode bound also suggests 
that r-mode spindown should be fast, i.e. at least compa- 
rable to magnetic spindown, which would require a high 
r-mode saturation amplitude a > 0.01 in young neutron 
stars. 

R-modes are an important source of gravitational 
waves and provide an interesting possibility for their de- 
tection since r-modes can radiate over long times. Pre- 
viously analytic estimates have been given for the grav- 
itational wave strain of r-modes of young stars [Q. As 
found in the latter work the most important quantity for 
the detectability of the gravitational wave signal is the 
minimum frequency emitted during the evolution. Here 
we extend this analysis and our semi-analytic results al- 
low us to determine the minimum gravitational strain 
expected from r-modes as well as the signal to noise ra- 
tio for their detection in forthcoming second generation 
detectors like advanced LIGO [9, 10 . Furthermore we 
point out that the limited observation time can strongly 
reduce the detectability of such sources. 
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II. MATERIAL AND STAR PROPERTIES 

As noted before a neutron star is a complex system 
and the description of its spindown evolution requires 
an understanding of various different aspects. The star- 
ing point is the static star configuration which is deter- 
mined by gravity as described by the general relativis- 
tic Tolman-Oppenheimer-Volkov (TOV) equations 
The latter require the equation of state of charge-neutral 
dense matter in beta equilibrium as an ingredient which 
is determined by microscopic strong interactions. Both 
the rotation and the oscillation of the star are modeled as 
perturbations of the static respectively the rotating star 
configuration |12| . For simplicity both are standardly 
obtained from the Newtonian Euler-equation of a non- 
viscous fluid. The classical r- modes with I = m are given 
by the Eulerian velocity perturbations 



5v = am 



(1) 



where a is a dimensionless amplitude parameter, R is 
the star radius, and lo are rotational and mode angu- 
lar velocity, respectively, and is the magnetic vector 
spherical harmonic with magnetic quantum number m. 

To describe the dynamical evolution requires micro- 
scopic material properties of the dense matter within the 
star. Here we follow the method that was developed pre- 
viously in and to obtain semi-analytic results for 
the boundary of the r-mode instability region as well as 
for the saturation amplitude in case suprathermal bulk 
viscosity saturates the mode. The basic idea is that the 
microscopic material properties that are relevant for the 
pulsar evolution have - at least over the range required to 
describe the relevant aspects - simple power law depen- 
dencies on the temperature T, the oscillation frequency 
Lo and the density oscillation amplitude An/h. The re- 
quired quantities for the star evolution are the shear vis- 
cosity rj, the bulk viscosity C, the specific heat cy and 
the neutrino emissivity e. Both the the bulk viscosity 
and the neutrino emissivity are determined by a weak 
equilibration process whose rate in a degenerate system 
has the form F^'*^' = —TT^fiA, where /lja is the chemi- 
cal potential difference that is driven out of equilibrium. 
Here we neglect the suprathermal amplitude dependence 
of the bulk viscosity and the neutrino emissivity [131 114] 
and leave this for future work. In the subthermal regime 
these quantities can generally be parameterized as 



cy, e, • • • . The important point is that the exponents 
are the same for a given phase of dense matter, like e.g. 
hadronic matter with modified Urea reactions, and quan- 
titative details due to the unknown equation of state or 
the interactions are only reflected in the prefactors. Strik- 
ingly we will see below that the quantitative dependence 
on these prefactors can be even far more insensitive than 
this general argument suggests. Different phases of dense 
matter, however, can feature qualitatively different low 
energy degrees of freedom resulting in different exponents 
and thereby can lead to a qualitatively very different star 
evolution. 

The above material properties are local quantities, but 
all that enters the evolution equations below are quanti- 
ties that are averaged over the entire star. These are the 
power radiated in gravitational waves Pg, the dissipated 
power Ps and Pb due to shear and bulk viscosity, the 
specific heat of the star Cv, the total neutrino luminos- 
ity and the momentum of inertia / O |S] 
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where Aqcd and Aew ^tre characteristic strong and elec- 
troweak scales introduced to make these quantities di- 
mensionless. With the energy of the r-mode 
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one can define characteristic time scales for the r-mode 
growth tg and the viscous damping ts and tb via 
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where B and C are non-equilibrium susceptibilities. 
With such power law behavior it is clear that physical 
results should depend far more sensitively on the expo- 
nents a, S, V and 6 than on the pref actor functions fj, 



so Tg is negative and ts and tb are positive. The above 
quantities are given in terms of parameters that involve 
the integration over the star or alternatively over the 
layer ranging from to R^ that contributes dominantly 
to this quantity. These are defined in table |T] 

All quantities appearing in these integrals, like the en- 
ergy density p, the inverse squared speed of sound A, the 
non-equilibrium susceptibility C and the r-mode density 
oscillation 
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Table I: Radial integral parameters encoding the complete 
information on the star's interior - i.e. on the microphysical 
transport properties, the equation of state and the star's den- 
sity profile. 
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depend on the position within the star via their density 
dependence. These few constants encode the complete 
information on the physics inside the star and are suffi- 
cient to calculate the star evolution. 

The dimensionless parameters / and J, arising in the 
moment of inertia and the canonical energy and angular 
momentum of the mode, involve integrals over the energy 
density in which the outer parts are strongly weighted by 
high powers of r. They are normalized by the mass and 
appropriate powers of the radius so that the range these 
constants can vary over is independent of the mass. The 
mass by which these constants are divided is given by 
an analogous integral with a smaller power of r. Conse- 
quently these constants are large for stars whose mass is 
strongly spread out but small for stars where the mass is 
concentrated close to the center. Since the energy den- 
sity is a monotonously decreasing function, it is clear 
that / and J are bounded from above by the values for 
a constant density star |3j. Similarly, lower bounds can 
be obtained by noting that in a star the mass cannot be 
arbitrarily concentrated. A limit for the energy density 
at a given radius is obtained by the constraint that the 
matter within this radius has to be stable against gravita- 
tional collapse. The corresponding bound for the energy 
density is given by p{r) < l/lSnGr'^) which can be inte- 
grated to obtain lower bounds on / and J. Combining 
these rigorous limits we find that for any compact star 
these parameters are narrowly bounded within roughly a 
factor of two 
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3.41 • 10-^ . (14) 
However, for a neutron star with a crust the energy 
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density vanishes at the surface and therefore the up- 
per bounds should even clearly overestimate the realistic 
range. 

The uncertainty on the other parameters arises both 
from the microscopic quantities eqs. ([2| and ^ as well as 
from the particular (baryon) density profile rijr) which in 
turn depends both on the equation of state and the par- 
ticular solution of the TOV equations (parameterized e.g. 
by the star's mass). All parameters defined above in tab. 
H] are given for different stars in tab. |ll] For illustration 
of our semi-analytic results below we consider different 
stars with an APR equation of state jT5]. In these stars 
we assume that the core of the star dominates the rel- 
evant quantities and neglect possible contributions from 
the crust. For the maximum mass APR star the density 
is high enough that direct Urea processes are kinemati- 
cally allowed within an inner core. 



III. PULSAR EVOLUTION 

The evolution equations are obtained from energy and 
angular momentum conservation laws [SJ[T5] and take the 
form 
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with Q = 3J/(2/) and in terms of the viscous dissi- 
pated power Pv = Ps + Pb + • • • and damping time 
1 /tv = 1/ts + 1/tb + ■ ■ ■ , where the dots denote pos- 
sible other dissipative mechanisms, like boundary layer 
effects. Here we include the reheating due to the dis- 
sipated power Pv within the star, that was partly ne- 
glected in [S]. Using the previous bounds on the param- 
eters / and J we see that Q < 81/(1127r) w 0.23 so that 
the factors 1 ± Qa^ are only relevant for large amplitude 
modes with a > 1, i.e. in a regime where the perturba- 
tive approximation for the r-mode |12| breaks down any- 
way. The r-mode is unstable when the right hand side 
of the amplitude equation eq. ( |15| is positive, i.e. when 
Tyj^^Q > —tg, where the equality defines the bound- 
ary of the instability region [3]. Since the r-mode will 
grow once it enters the instability region the evolution re- 
quires a non-linear mechanism to saturate the amplitude 
at a finite value. A simple possibility is the suprathermal 
enhancement of the bulk viscosity [H [131 Hi] • Unfortu- 
nately when only the damping in the core is considered, 
this mechanism saturates the mode only at rather large 
amplitudes. The bulk viscosity contribution of the in- 
ner crust might change this since the outer regions of the 
star are strongly weighted by the suprathermal viscosity 
[U , but the bulk viscosity has not been computed at this 
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Table II: Parameters characterizing the neutron star considered in this work. The constants I, J, S, V, Cv and L are given in 
tab. |l] using the generic normalization scales Aqcd ~ 1 GeV and Aew = 100 GeV and the temperature exponents a, S, v and 
6 are defined by eqs. (|2| and ([sjl. 



point. Another promising mechanism is the non-linear 
coupling of the r-mode to other daughter modes that are 
subsequently damped by viscosity [I84. Further possibil- 
ities include non-linear hydrodynamic effects |19[ [20] or 
turbulent crust-core boundary layer damping |21| but at 
this point it is not settled which of these mechanisms will 
actually saturate the mode. Therefore, we follow here the 
approach pursued in j8| and will not study a particular 
saturation mechanism but simply assume that a corre- 
sponding mechanism operates and saturates the r-mode 
at a particular amplitude asat- In contrast to the above 
explicit mechanisms, asat is in this case an unknown pa- 
rameter that may depend on 51 or T and we will study 
the dependence of our results on it below. 
At saturation da/dt = from eq. (15), 
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and likewise Py — > Pq. This simply says that at satu- 
ration there must be a sufficiently strong source of dissi- 
pation to overcome the gravitational instability. Making 
the replacements in eqs. (151, (16 1 and (17) leaves the 
reduced set 



can be a general function agat {T, This function has 
to vanish at the boundary of the instability region. Due 
to the strong power law dependence of the corresponding 
time scales tq and Ty on T and fJ, it drops sharply over 
a very narrow region, as has been explicitly shown in [4J. 
For the spindown evolution this narrow boundary region 
is far less important than the much larger interior of the 
instability region. In the interior the amplitude will de- 
pend much more weakly on T and J7 and we make the 
simple power law ansatz 



(21) 



For instance the semi-analytic expression for the satu- 
ration amplitude due to suprathermal bulk viscosity has 
such a power law dependence on the frequency and proves 
to be even independent of temperature, i.e. /? — [4J. 
Since the only temperature dependence in eq. ( 19 1 comes 
from the implicit dependence via a^ati for /3 = the spin- 
down evolution is not affected at all by the thermal evo- 
lution. This applies in particular to the model employed 
in |8] where the saturation amplitude is assumed to be 
constant throughout the instability region. 

To analyze the above equations it is useful to introduce 
characteristic evolution time scales 
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where, as argued before, the approximate expressions 
hold over nearly the entire range of physically reasonable 
amplitudes. Note that compared to [8 we take into ac- 
count that the saturation mechanism dissipates according 
to eq. ( 18 1 a significant amount of energy that heats the 
star. The only way that this could be avoided is if some of 
the corresponding energy is radiated away, e.g. in neutri- 
nos or electromagnetic radiation. In all proposed mecha- 
nisms, like the non-linear enhancement of the bulk viscos- 
ity lU [T7] , mode-coupling to viscously damped daughter 
modes [TB] or nonlinear hydrodynamic effects |19l 120] . 
this is not the case and the dissipated energy eventually 
entirely ends up as heat. 

The complete information on the saturation mecha- 
nism is then encoded in the saturation amplitude which 



(22) 

A young compact star is born in a core bounce super- 
nova with a large temperature T > 10^^ K. We will 
assume that the initial frequency is above the mini- 
mum frequency of the instability region [3 . Initially 
Tx ^ j^-{o-v-i) jg sjiiall so that the star cools very fast 
due to neutrino emission and quickly enters the insta- 
bility region. There by definition [tq[ < tv[^^q and 
(while the star cools further) the r-mode will grow with 
time scale Tq, « tq. For a compact star spinning with 
a millisecond period this scale is of the order of sec- 
onds, so that the amplitude quickly grows until it is sat- 
urated at some Ofgat- From this point on the fast evolu- 
tion of a stops and it changes slowly with temperature 
and frequency. Since cooling is initially very quick but 
slows down considerably at lower temperatures, at some 
point the temperature-independent dissipation eq. ( 18 1 
required to saturate the mode will dominate. The im- 
portant question for the evolution is then if the spin- 
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down or the temperature change is faster. To see this let 
us analyze the behavior when temperature change due to 
dissipative heating (neglecting neutrino emission) Tj, is 

faster than the spindown, i.e. < tq. This condi- 

tion yields a relation for the angular velocity of the star 
which introducing the Kepler frequency fl^ = |'\/27rG'/3o 
can be written in the form 




v + l 



> const-— 



(23) 



where in the second line we use the fact that for dense 
matter in general the temperature exponent arising in 
the specific heat eq. (|3| obeys v > 1, and larger val- 
ues make the right hand side even smaller. For a com- 
pact star the constant is 0(1), since the dimensionless 
parameters Cy and / are normalized to also be 0(1), 
the characteristic scale for the average energy density 
p of strongly interacting matter is Aqqj-,, the radius 
of a neutron star is close to its Schwarzschild radius 
Rs = 2GM and the root further decreases deviations 
from unity. However, since a compact star is a degener- 
ate system, after a short time it reaches a temperature 
where T/Aqcd ^ O{l0^'^). Therefore, taking into ac- 
count the O (1) prefactor, the viscous heating should be 
faster than the spindown for frequencies il/flK ^ 10~^. 
Since r-modes are only present within the instability re- 
gion, this should be compared to the minimum of the 
instability region Omin, for which a semi-analytic expres- 
sion is given in [I?. There it was shown that this mini- 
mum is extremely insensitive to the microscopic details 
which yields a general limit of roughly ft^^ij^/ftx 1/20. 
Therefore, at r-mode saturation within the instability re- 
gion the dissipative heating is always faster than the spin- 
down. Note that this conclusion is generic and to leading 
order independent of both the saturation amplitude and 
the composition of the star. 



The evolution should therefore follow this curve until 
the star leaves the instability region. As can be seen in 
fig. [T] the boundary of the instability region of a neutron 
star has different segments: a low-temperature bound- 
ary outside of which shear viscosity damps the mode and 
an intermediate-temperature boundary where bulk vis- 
cosity is the dominant damping mechanism. Since the 
bulk viscosity features a resonant behavior and decreases 
again at very large temperatures T > 10^^ K, there a 
is third high-temperature boundary above which the r- 
mode is unstable (see fig. 3 of [3]), but for neutron stars 
this segment is quite likely physically irrelevant since the 
stars cools below it quickly. Semi-analytic expressions 
for these segments of the boundary of the instability re- 
gion have been derived in [3 . The form of the boundary 
depends strongly on the form of matter. When exotic 
forms of matter are present then in the relevant temper- 
ature regime there are so called stability windows, caused 
by the resonant behavior of the bulk viscosity of the cor- 
responding phase, where the star is stable against devel- 
oping an r-mode up to large frequencies (see fig. 3 of [3]). 
In case the steady-state curve formally lies in such a sta- 
bility window it would clearly be irrelevant for the spin- 
down^. Then the evolution can take a long time, since 
the strong heating can push the star out of the instability 
region before the r-mode reaches a large amplitude, and 
other spindown sources will likely dominate. Therefore 
we will restrict ourselves here to the case of pure neutron 
stars. The general evolution for different forms of matter 
will be discussed elsewhere. 

In the case of neutron stars the intermediate- 
temperature boundary of the instability region has a 
steeper slope than the steady state curve which thereby 
can only intersect the low-temperature boundary. In 
the case of a neutron star the low-temperature bound- 
ary Qib(T) is determined by the condition ts = tq and 
the corresponding semi-analytic solution is [5] 
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IV. SEMI-ANALYTIC EXPRESSIONS 

A. Endpoint of the evolution 

Since for a compact star the thermal evolution is faster 
than the spindown, the evolution of the star should follow 
a curve ^hc{T) where it is always thermally in a steady 
state where neutrino cooling equals viscous heating + 
Pq = fTHl. With the semi-analytic expressions eqs. Q 
and ([s]) this equation can be solved explicitly and yields 
for the dominant m = 2 r-mode 



The final angular velocity f2 j where the evolution leaves 
the instability region is given by 
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^ As just noted, for standard neutron stars the corresponding sta- 
bility window lies at very large temperatures T > 10^" K where 
the neutrino cooling still strongly dominates. 
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and the corresponding temperature Tf reads 



Tf = 




2+27 
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Note that, if we had used a different steady-state curve 
i s) 

^hc where the heating conies only from shear viscosity, 
as assumed in [8J , or if we had included the bulk viscosity 
as well, the rest of the evolution could be very different, 
but we would have obtained exactly the same result for 
the endpoint eqs. (26 1 and (27 1 of the evolution . This 



surprising feature comes about since the shear viscosity 
is strongly dominant compared to the bulk viscosity on 
the left boundary and this boundary is precisely defined 
by the condition T5 = tq which implies that dissipated 
power equals the power radiated in gravitational waves 
Ps^Pq- In this sense the result eq. (26 1 is universal and 



does not depend on assumptions about the dissipative 
aspects of the saturation mechanism. 

The above expressions involve the star constants given 
in tab. [IT] In order to give these expressions in a more 
physical form we write them relative to a fiducial refer- 
ence model chosen as a 1.4 neutron star with an APR 
equation of state [TS] given in the first row of tab. [ll|and 
denoted by a subscript "fid". In case of a constant, T- 
and f2-independent, saturation amplitude [8 f3 — j — 
this gives for a standard neutron star (NS) with shear 
viscosity due to leptonic processes and neutrino emission 
due to modified Urea reactions 
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First it is interesting to note that all these expressions 
depend neither on the specific heat nor directly on the 
bulk viscosity of the considered form of matter. The 
latter is at first sight surprising since the bulk viscosity 
is the dominant dissipation mechanism at high temper- 
atures and low amplitudes. But during the spindown 
evolution the dissipation is dominated by the appropri- 
ate saturation mechanism while the amount of dissipa- 
tion is according to eqs. (11) and (18 1 determined en- 



tirely by gravitational physics. This is particularly fa- 
vorable since it eliminates the complications arising from 



the fact that the damping due to bulk viscosity requires 
the density fiuctuation of the r-mode eq. ( 12 1 which is 



only non- vanishing to second order in Q. [T5] and involves 
significant uncertainties. Even more striking is the ex- 
treme insensitivity of these result to the remaining mi- 
crophysics included in the constants S and L. The neu- 
trino emissivity constant L, in particular, depends on 
poorly known strong interaction corrections to the weak 
processes Yet, the remarkably small power 5/184 

arising in the important expression for the final angular 
velocity fi^^'^-* almost eliminates this uncertainty. The 

dependence of fJ^^'^'' on the shear viscosity constant S 

is somewhat stronger but S is dominated by electromag- 
netic leptonic processes [53] that are well under theoret- 
ical control^ and depends on the equation of state only 
via the lepton densities. The dependence on the param- 
eters J, M and R which encode macroscopic properties 
of the star is stronger, but they vary only within narrow 
margins. This insensitivity to the underlying parameters 
is similar to the case of the minimum of the instabil- 
ity region analyzed previously in [3], yet it is fortunately 
even more pronounced in the physically important case 
eq. (29 1. The most striking feature of these expression, 
however, is that the insensitivity extends even to the 
dependence on ofsat- which is theoretically uncertain by 
many orders of magnitude [?1 ll8H20j and therefore rep- 
resents the main uncertainty in the analysis. Here the 
similarly small power 5/92 dramatically reduces this de- 
pendence and thereby allows a quantitative comparison 
with observational data below. 

To get an idea of the uncertainty in fj^^^ we esti- 
mate the error ranges of the individual parameters by 
symmetric logarithmic error bands - i.e. we allow for the 
corresponding parameter to be smaller or larger than the 
APR reference star by a given factor. The shear viscosity 
parameter S due to leptonic processes [53] is uncertain 
within a factor of two of Sud and the neutrino emissivity 
parameter L due to modified Urea processes |22| within 
an order of magnitude of ifid- For J we have rigorous lim- 
its eq. ( 14 1 . An upper limit for the mass of a neutron star 



is < 2.5 Mq and an upper limit for the radius of a star 
rotating with millisecond periods about < 15 km. Since 
the upper bound for J is for a constant density profile 
which is far from the situation in a neutron star and, due 
to cancellations, the dependence on the mass and radius 
turns out to be significantly weaker than these individual 
individual errors suggest, see table |IIl| below, we consider 
half the above ranges for the individual logarithmic un- 
certainties for J, M and i? as a more realistic estimate. 
Under the assumption that these errors are independent 
we find 



^ There is a hadronic contribution to the shear viscosity as well 
but the latter is strongly subleading at temperatures relevant to 
young neutron stars due to the different power law dependence. 
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(61.4 ±9.4) Hz a, 



-5/92 



(30) 



We want to stress that this is just an estimate of the 
Hkely error range and not a statistical measure of well 
defined significance. 

Different classes of compact stars, containing other 
phases of dense matter with different low energy degrees 
of freedom, feature different values for the exponents in 
eqs. ^ and (|3| and can lead to qualitatively different be- 
havior. An example is a heavy star where direct Urea re- 
actions are kinematically allowed in an inner core and the 
parameters change according to table |TT] which increases 
the final frequency beyond the above error range. When 
exotic forms of matter are present, like hyperon or quark 
matter, the enhanced dissipation of these phases leads to 
stability windows |31 [23] . In this case the boundary of 
the stability window in the relevant temperature regime 
is determined by the bulk viscosity of the exotic phase, 
so that eqs. (26 1 and (27 1 are not valid. The steady state 



curve then ends at a higher frequency due to the insta- 
bility window where the star leaves the unstable regime 
and cools until it reaches the lower boundary of the sta- 
bility window. Since at these lower temperatures r-mode 
heating dominates neutrino cooling the star is repeatedly 
pushed out of the instability region once its amplitude be- 
comes large and the residual spindown takes much longer 
[251 The important point is however, that the final 
frequency which is in this case given by the minimum 
of the instability region and for which a similarly pre- 
cise semi-analytic result exists |i3j is systematically larger 
for exotic forms of matter since the strong bulk viscosity 
of the exotic phase dominates the shear viscosity down 
to lower temperatures. Therefore, the above expressions 



eqs. (29 1 and (p50|) provide lower bounds for the frequency 



to which an arbitrary compact star can be spun down due 
to the r-mode instability. 



B. Evolution time scales 

Next let us actually solve the evolution equations to 
obtain the corresponding evolution time scales. The ini- 
tial cooling time in the stable regime before the evolution 
enters the instability region tg^ results from solving the 
thermal equation eq. (17 1, in the absence of damping, 

letermined by the semi-analytic 



down to the temperature d( 
expression for the intermediate-temperature boundary of 
the instability region given in [5]. The result reads 



1 



3c2 



3^5 



- w - 1 V 2i277rK2 



(31) 



1/5 



QCD 



-v-1 



4(e-i5— u-l) 
EW 



pGh/PVLf 



in terms of the initial angular velocity fi^ and depends 
both on the bulk viscosity and the specific heat. For a 
standard neutron star with shear viscosity due to leptonic 
processes and neutrino emission due to modified Urea 
reactions this gives 



.{NS) 



7.83-10 




This time is so short compared to the subsequent r-mode 
evolution that it can safely be neglected. 

The spindown time tgd is obtained by solving the fre- 
quency equation eq. (19) in the saturated regime. As 



noted before, this equation completely decouples for /3 = 
and factoring the frequency dependence out of the grav- 
itational time scale via tq = tg^~^ yields the frequency 
solution 



VL{t)= \ VL 



-6-27 



4(3 + 7) Q«L^ 



-1/(6+27) 



(33) 



We find in terms of the gravitational time scale eq. ^ 
evaluated at the minimum frequency 



l-sd 



12(1 



-7/3)g«L 
1 



1 - 



6+27 N 



nf^fii 6(l + 7/3)ri 



(34) 



/ 



If the saturation amplitude has a frequency dependence 
with exponent 7 < 0, then up to small corrections in 
flf /Cli this simply increases the spindown time com- 
pared to a constant saturation amplitude by a factor^ 
1/ (1 -|- 7/3). In case of dissipation due to suprathermal 
bulk viscosity (1] the spindown time is longer compared 
to the saturation model with a constant amplitude by a 
factor of 9/5. 



The above expression eq. ( 34 1 looks similar to the rela- 



tion between the spindown time for other spindown mech- 
anisms, like magnetic dipole radiation or gravitational 
wave emission due to deformations, and the character- 
istic pulsar time scale — Oo/i^o 13- Yet, the important 
point to note here is, that Qf and ilf are the values be- 
fore the star leaves the instability region which can be 
very different from the observed values JIq E^nd flo if the 
observation is performed after the star has left the in- 
stability region. The characteristic feature of the r-mode 



^ In the case 7 = —3 the spindown equation eq. \l9\ would give a 
logarithmic instead of a power law solution, yet this special case 
does not seem physically relevant. 



8 



mechanism, namely that it operates only for a finite time 
interval, renders the r-mode spindown time completely 
independent of the characteristic pulsar time scale and 
allows the former to be orders of magnitude smaller, since 
the spindown rate can strongly drop from a large r-mode 
rate to a much lower rate of other spindown mechanisms 
at the boundary of the instability region. 



Inserting the expression for the final frequency eq. ( 34 1 



into eq. ( 29 ) this yields in case of a standard neutron star 
and for ilf <C f^i 



18/23 



JNS) 




-15/92 



-107/92 




-5/46 



-77/46 



(35) 



Here the dependence on the parameters is stronger so 
that the duration is considerably more uncertain. In 
particular it is very sensitive to the saturation ampli- 
tude Ofsat- If Q^sat becomes too small this time scale is 
so long that other spindown mechanisms will dominate, 
as is clear from the sizable observed spindown rate of 



young pulsars [57]. Employing the bounds eqs. (13) and 



( 14 1 as well as the full uncertainty ranges for the other 



underlying parameters discussed before to estimate the 
uncertainty in the final frequency eq. ( 30 1 we find that 



the spindown time could be smaller or larger by nearly 
an order of magnitude. 



V. SPINDOWN OF YOUNG PULSARS 

A. Comparison of analytic and numeric results 

The results for the evolution of the 1.4 neutron 
star in temperature-angular velocity space is shown in 
fig. [T| for different values of the initial rotation angular 
velocity and the saturation amplitude. As can be seen 
the numerical solution given by the solid curves perfectly 
confirms the above picture. The star initially cools until 
it reaches the instability region, where the r-mode ampli- 
tude starts to grow. Once the amplitude is significant 
the star starts spinning down but since the spindown 
time tq is much longer than the cooling time tt, the 
star cools further until it reaches the steady state curve** 
where the fast temperature evolution is stopped and the 
star evolves along the steady state curve on the longer 
time scale tq. Consequently, the evolution is completely 



* For the largest saturation amplitudes the evolution reaches the 
steady state temperature before its amplitude has saturated. 
Therefore, as can be seen from the Osat = 1 curve in fig. [l] the 
star "undercools" to lower temperatures and then reheats back to 
the steady state curve once the saturation amplitude is reached. 




Figure 1: The spindown evolution of a young 1.4 Mq neu- 
tron star with an APR equation of state |15| in temperature- 
angular velocity space. The boundary of the instability region 
of the fundamental m = 2 r-mode is shown by the dotted 
curve. The dashed curves (which are mostly hidden under- 
neath the solid curves) represent the steady state eq. ( |24[ ) 
where heating equals cooling and are given for different r- 
mode amplitudes ranging from Qsat ~ 10"'' (left) to dsat = 1 
(right). The solid lines show the numerical solution of the 
evolution equations for two fiducial initial spin frequencies 
n — 0.8 ilx and Q = 0.2 ^1k- As can be seen, after an initial 
cooling phase the evolution simply merges on the appropriate 
steady state curve and follows it to the edge of the instabil- 
ity region. The dots denote the semi-analytic results for the 
points of the spindown evolution wh ere t he s tar enters and 
leaves the instability region, eqs. ( 28 1 and ( 29 1 . 



independent of the initial amplitude amin required to 
start the growth in eq. (15 I. Such initial perturbations 



will clearly be present immediately after the creation of 
the neutron star following a core bounce. As anticipated, 
the numerical solution also proves that the evolution is to 
very good accuracy independent of the initial frequency 
unless it is very close to the final frequency. Therefore, 
similar to the thermal evolution, the spindown evolution 
completely loses its memory of the initial conditions at 
sufficiently late times. The saturation amplitude in con- 
trast provides the dominant uncertainty for the evolu- 
tion and the trajectory is shown for a range of possible 
saturation amplitudes. Larger values lead to stronger 
dissipative reheating and therefore spin down the star 
at larger temperatures so that the evolution leaves the 
instability region at lower frequencies closer to the min- 
imum of the instability region. The dots are obtained 
from the analytic expressions eqs. (28 1 and (29 1 which 
agree with the numerical solution. Note again that we 
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neglect other spindown mechanisms Hke magnetic dipole 
radiation in this work which would speed up the evolu- 
tion and continue the magnetic spindown after the star 
leaves the r-mode instability region. 

In jE] the evolution was artificially stopped at a temper- 
ature of 10^ K since superfluidity and non-perfect fluid 
effects were expected to invalidate the analysis below. 
We do not follow this ad hoc prescription here since it 
is not likely that these effects have such a big influence 
on the spindown evolution for the following reasons: One 
important effect of superfluidity is to affect beta equili- 
bration. Around the critical temperature it is enhanced 
due to pair breaking, but at temperatures far below and 
small oscillation amplitudes beta equilibration processes 
are strongly suppressed. Since the superfluid gap de- 
pends strongly on density the entire star will not be 
gapped at such temperatures and therefore in part of the 
star the weak interactions are hardly affected and the 
corresponding parameters L and V are only moderately 
reduced. Even more, as has recently been shown in |14| . 
superfluidity does not suppress weak processes at all at 
sufficiently large amplitudes that are within the range we 
study here. Superfluidity can also lead to an additional 
source of dissipation via mutual friction. In contrast to 
the earlier generic analysis referred to in [8J, the later 
dedicated r-mode analysis [2S] showed that it is unlikely 
that mutual friction has such a dramatic impact on the 
r-mode instability. In |.29j this effect was studied, taking 
into account the large uncertainties in pairing and the mi- 
croscopic drag parameter. It was found that depending 
on these uncertainties the effect can range from the un- 
likely extreme that the r-mode is completely suppressed 
below a certain temperature to a nearly negligible impact 
on the instability region compared to the ungapped case. 
Therefore, it is interesting to study the present case of 
standard dissipation sources in detail, since it sets the 
limit on the frequency to which r-modes can spin down a 
star as well as on the emitted gravitational wave signal, as 
discussed below. Any additional dissipation source will 
reduce the r-mode instability and correspondingly lead to 
larger flnal spin frequencies and a smaller gravitational 
wave signal. 

Fig. |2] shows the evolution of the angular velocity as 
a function of time. As predicted by eq. (351, the spin- 



1.0 



down time increases strongly with decreasing saturation 
amplitude. The points show that the semi-analytic re- 
sults for this time again agree with the numerical solu- 
tion^. Here we are mainly interested in the spindown 



^ Note that he spindown overshoots at large amplitude since it 
takes time for the amplitude to decay and the spindown continues 
during this period, as can be seen from the lowest curve which 
drops below the dot so that the final frequency lies somewhat 
below the instability region. This is mainly an artifact of the 
simplified model which keeps the amplitude constant until the 
instability boundary is reached. In reality with an explicit satu- 
ration mechanism the amplitude has to vanish at the boundary in 
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Figure 2: The spindown evolution of a young neutron star. 
The solid lines show the evolution of a 1.4 M© star with an 
APR equation of state and for different amplitudes ranging 
from Osat ~ 1 (left) to asat ~ 10"^ (right). The dots denote 
the analytic results eq. ( 32 1 respectively eqs. ( 35 1 and ( 29 1 for 



the point of the spindown evolution where the star enters and 
leaves the instability region. 



aspects and only note that as far as the thermal evolu- 
tion is concerned the r-mode effectively causes a delay 
in the cooling. This is signiflcant initially, but, since the 
cooling slows down strongly, it is irrelevant at later times 
t > t,d- 

Table |III| lists the semi-analytic results for a range of 
saturation amplitudes for two different neutron stars - 
the previously discussed standard star with a canonical 
mass of 1.4 Mq and, motivated by the recent observation 
|30| . also a heavy 2.0 Mq star. Strikingly the flnal fre- 
quency proves nearly independent of the mass, and is far 
less sensitive than eq. (29) naively suggests. This is ap- 



parently caused by cancellations due to the implicit mass 
dependence of the radius and of the parameters J, S and 
L. Despite the explicit normalizations in table|l]these pa- 
rameters still have a signiflcant implicit dependence on 
M and R via the density proflle of the star. Correspond- 
ingly, whereas the dependence on microscopic material 
properties is explicitly given by the semi-analytic expres- 
sions eqs. (26 1 to (35 1, the dependence on star properties 



like the mass and radius is not properly reflected by the 
explicit dependencies alone but the additional implicit 



the static limit, so it must start decreasing inside the instability 
region, see |4]. 
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Q-sat 


Tf [10»K] 




// m 


isd [y] 


1 


12.6 (13.9) 


0.064 (0.049) 


61.4 (60.3) 


12.3 (9.4) 


0.1 


8.02 (8.88) 


0.073 (0.056) 


69.5 (68.4) 


582 (445) 


10"^ 


5.11 (5.66) 


0.082 (0.063) 


78.8 (77.5) 


2.75 (2.10) • 10" 


10"^ 


3.26 (3.61) 


0.093 (0.072) 


89.3 (87.8) 


1.30 (0.99) • 10*^ 


10-* 


2.08 (2.30) 


0.106 (0.081) 


101.2 (99.5) 


6.11 (4.68) • lO'^ 



Table III: Dependence of the results obtained from the semi- 



analytic expressions for the characteristic quantities eqs. ( 28 1 



( 29 1 and (351 of the spindown evolution for neutron star with 
an APR equation of state on the saturation amplitude Qsat. 
We give in each case results for both a standard 1.4 M0 star 
and a heavy 2.0 Mq star next to it in parentheses. At the 
smallest amplitudes the given spindown times tsd should not 
be relevant since they would be so long that the neglected 
magnetic spindown should dominate and spin down the star 
on much shorter time scales. 



dependence via the averaged quantities in table |T] is cru- 
cial. The initial cooling time before the star enters the 
instability region is amplitude independent and yields 
for a star spinning at the Kepler frequency the small 
values tsc = 96.4 (55.1) ms, which rise to a few hours 
for frequencies at the lower end of the instability region, 
which is indeed negligible compared to the long spindown 
times. The left part of table |IV| finally compares the spin- 
down observables for different stars at a fixed amplitude 
Q!sat = 1- We see that, when direct Urea processes are al- 
lowed, the final temperature is significantly lower due to 
the enhanced neutrino cooling. Since the instability re- 
gion shrinks at lower temperatures the final frequency is 
higher and the spindown time is correspondingly shorter. 



B. Saturation model dependence 



Let us now discuss the influence of different satura- 



tion models asat(r, Vl). Since in the general case eq. (21 ) 
the amplitude varies over the instability region and the 
coefficient dsat is a dimensionful parameter, we have to 
find analogous parameter values for dsat to compare sat- 
uration models with different values of /3 and 7 among 
each other and with the constant model discussed so far. 
Since the evolution spends most of the time in the final 
stages of the spindown evolution it is natural to choose 
Q^sat so that different models have a similar amplitude in 
this stage. If we compare models defined by dsat, P and 7 
that have the same amplitude Ofgat at the end of the evo- 
lution, then they leave the instability region at the same 
point given by eqs. (26 1 and (27 1. This can be easily 
seen, noting that the conditions T5 = tq and Ps — Pg 
that determine the endpoint depend only on the satu- 
ration amplitude aaat{Tf ^^f) at this point. However, 
the path within the instability region, determined by the 
steady state curve eq. ( p4| , and in particular the thermal 
behavior will depend on 13 and 7. 

Fig. [3] shows a few examples of different saturation 



models for two different saturation amplitudes asat = 
IQ-^ and asat = 1 at the end of the evolution. The 
solid lines are the constant model /3 = 7 = |H] dis- 
cussed above. The dashed curves show the case when 
suprathermal viscosity is strong enough to saturate the 
mode. In this case the analytic result for the amplitude 
is ckgat ~ Vr'^/^ and the amplitude becomes large at 
low frequency where the suprathermal enhancement of 
the viscosity is weaker. Since the amplitude is lower at 
large frequency the dissipative heating is smaller and the 
steady state curve moves to lower temperatures. Since 
the spindown evolution is in this case still independent of 
temperature, the spindown time eq. ( |35[ ) does not change 
compared to the constant model and the only difference 
is in the thermal evolution within the instability region. 
The thermal evolution could further be altered by the 
suprathermal enhancement of the neutrino emissivity |14| 
that we neglect here. 

The final example is a case where the amplitude rises 
with temperature asat ^ T . This seems to be realized to 
some extent in mode-coupling models, where the satura- 
tion amplitude is as low as asat — O (lO"^) in the low 
temperature regime relevant for low mass X-ray binaries 
[31 which is smaller than in the high temperature regime 
relevant for young stars where asat = O(l0-^ — lO"^) 
|18] . As shown by the dot-dashed curves {/3 = 1, 7 = 0) 
in fig. [3] the steady-state curve lies at higher tempera- 
tures^. Because of the temperature dependence of the 



saturation amplitude in eq. ( 19 1 this could now actually 



change the spindown time. However, there is only a sig- 
nificant increase of the amplitude at high temperatures 
and thereby higher frequencies. Since the evolution is 
fast in this region anyway this would hardly affect the 
overall spindown time. Therefore, we do not expect a 
qualitative change of the picture presented above in this 
case either^. 

In conclusion, different saturation mechanisms could 
change the thermal evolution but should not change the 
spindown behavior within the uncertainty range for the 
saturation amplitude from the picture obtained above 
with a constant amplitude. Furthermore, eq. (26 I shows 
that the observed insensitivity to the microscopic param- 
eters is not strongly altered by the saturation model for 



In mode coupling models the saturation amplitude also has a fre- 
quency dependence with a negative exponent which would coun- 
teract the temperature effect and bring the curve closer to the 
constant model again. 

At large saturation amplitudes the steady state curve crosses 
the instability boundary at an angular velocity below the Kepler 
frequency (see the rightmost curve in fig. |3]l. This could ni 
principle result in a qualitatively different initial evolution that 
oscillates around the boundary since strong reheating is pushing 
the evolution out of the instability region. Yet in this case the 
amplitude at high temperatures is unphysically large. Note that 
the amplitude is limited by energy constraints to values < O (1) 
since the mode energy cannot become bigger than the rotational 
energy. 
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spindown models. The time derivative of the frequency 
is readily obtained from eq. (19). Fig. |4] compares the 
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Figure 3: The dependence of the steady-state curve eq. (24 1 
on the saturation model for two different saturation ampli- 
tudes Osat {Tf,Q,f) — 10~* (left curves) and 1 (right curves) 
at the end of the evolution. The solid line shows the model 
with constant amplitude /3 = 7 = [8], the dashed line de- 
notes the case = 0, 7 = —4/3 corresponding to saturation 
by non-linear bulk viscosity [4. and the dotted curve shows a 
fiducial model with /3 = 1, 7 = which might capture certain 
features of the saturation model via mode coupling |18l I31| . 



realistic values 



I7I < 



C. Comparison to pulsar data 

The semi-analytic expressions finally allow us to com- 
pare our theoretical results to observational data. In 
most analyses [32J [33], merely the observed pulsar fre- 
quencies are compared to numerical computations of the 
r-mode instability region. For stars which spin slow 
enough that they are clearly below the minimum of the 
instability region this shows that they cannot have r- 
mode oscillations at present. However, there is at least 
one young star which lies above the minimum of the insta- 
bility region. Since the instability boundary is strongly 
temperature dependent, information on the temperature 
of fast stars is required to decide if they are actually inside 
the instability region. Yet, the temperature is unknown 
for most pulsars. Therefore, we will here go beyond such 
a static analysis and take into account the dynamic evo- 



lution as well as the frequency eq. (26 1 where the star 
actually leaves the instability region. 

In addition to the spin frequencies also the spindown 
rates are known for many young stars. Naturally, this 
additional information provides a stronger constraint on 



spindown evolution for different saturation amplitudes 
to the data for young stars that still have a significant 
spindown rate. As can be seen there is one pulsar J0537- 
6910 that spins fast enough that it could be within the 
instability region of a standard neutron star whose min- 
imum frequency is denoted by the vertical dotted line to 
the left^. Taking into account the dynamical evolution, 
the additional knowledge of the spindown rate alleviates 
the uncertainty with respect to the saturation amplitude, 
since larger amplitudes lead to significantly larger spin- 
down rates. This is shown by the black line which is 
formed by the semi-analytic results for the endpoints of 
the evolution and which represents the boundary of the 
instability region in /-/-space. The gray band reflects 
the uncertainty in the underlying parameters eq. (30 1. 



The fastest pulsar J0537-6910 is obviously compatible 
with being outside of the instability region. Thereby no 
observed star lies definitely inside the instability region 
which confirms the r-mode scenario for a sufficiently short 
spindown time. Moreover, J0537-6910 lies clearly below 
the line for our APR neutron star and just at the lower 
edge of the uncertainty band. Therefore, very likely no 
young star observed so far currently emits gravitational 
waves because of the r-mode instability. This is strength- 
ened by the absence of gravitational wave signals from 
observed young pulsars - including in particular J0537- 
6910 - in the latest LIGO data [34]. 

Our analysis was performed for a standard neutron 
star, but the instability regions for other possible exotic 
forms of matter do not extend to such low frequencies 
and thereby cannot change the conclusion that the ob- 
served stars are out of the instability region. Due to the 
insensitivity of the final frequency eq. ( p6| to the un- 
derlying parameters this should be a robust result and 
could be explicitly confirmed once a braking index for 
J0537-6910 can be measured. As a consequence our re- 
sult shows that the r-mode picture of the spindown of 
young pulsars presents a viable explanation for the low 
spin rates of young pulsars if the saturation amplitude 
is large enough ctsat > 0.01 to spin down the pulsar in a 
time tg^ < 10* y. An upper limit for the saturation am- 
plitude c^sat < 1 comes from the observed timing data of 
J0537-6910, as noted in [3S], since r-modes would have 
otherwise spun it down to its current frequency in less 
than a hundred years which is much shorter than its age 
estimate [35] . Whereas the r-mode picture yields a quan- 
titative explanation for the maximum observed spin fre- 
quency, other mechanisms, like the possibility that all 
stars are already born with lower frequencies or the sub- 
sequent spindown via magnetic dipole radiation, have as 



* The analytic result for the minimum of the instability region 
given in [3] would allow us to estimate also the uncertainties 
on the minimum frequency, but we refrain from this for better 
readability. 
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Figure 4: The spindown evolution of a young neutron star 
compared to observed pulsar data from the ATNF catalog 
|27| . The solid lines show the evolution of the 1.4 Mq star for 
different saturation amplitudes ranging from ctsat = 1 (top) 
to Qfsat = 10~* (bottom). The vertical dashed line to the 
left shows the lower frequency boundary of the instability re- 
gion. The steep solid line is formed by the endpoints of the 
evolution for different values of a and thereby represents the 
boundary of the instability region in /-/-space, whereas the 
gray band reflects the error due to the uncertainty in the 
underlying parameters. The dotted line segments show the 
current evolution for two stars for which a reliable braking 
index is available [5] • 



far as we know not provided a robust explanation why 
fast stars with a small initial magnetic dipole moment 
resulting in a slow spindown could not form. 



VI. GRAVITATIONAL WAVE EMISSION 

As discussed in the previous section, all presently ob- 
served young pulsars are likely outside of the instability 
region and so should not emit gravitational radiation due 
to r-modes. Nevertheless, it is possible that in the future 
we may observe a young pulsar that is still in the unsta- 
ble region. Furthermore, gravitational waves might even 
be detected from so far unknown young compact stars 
that have not been seen because their electromagnetic 
radiation is either too faint or it is absorbed or outshined 
by the supernova remnant. This possibility is interesting 
since there are several known young supernova remnants, 
including most notably SN 1987A, where a compact ob- 
ject has not been observed, yet. In view of the upcoming 
next generation of gravitational wave detectors like ad- 
vanced LIGO [9 it is therefore important to understand 



the expected gravitational radiation of such sources in 
detail. 

Above we already gave semi-analytic expressions for 
two important quantities in this respect, namely the ro- 
tational angular velocity il/ eq. (26 1 when the star ex- 
its the r-mode instability region, which is related to the 
lower bound Vf of the gravitational wave spectrum via 
V = 21 (Stt) 57, and the spindown time tsd eq. (34 1, which 



is the maximum age up to which gravitational waves from 
such a source could be observed. For a standard neutron 
star eq. (30 1 we estimate Vj « (81.9 ± 12.5) Hz a. 



-5/92 



To have a reliable estimate for the lower bound of the 
spectrum is particularly important since advanced LIGO 
will have a significantly improved sensitivity below 100 
Hz which coincides with the final stages of the r-mode 
spindown. Since the r-mode evolution becomes slow at 
low frequencies this region is where the evolution spends 
the most time and the analytic expression for the spin- 
down time function of the dependent parameters 
should then allow us to estimate the probability that 
sources within the reach of advanced LIGO are currently 
active and observable. 



A. Gravitational emission due to r-modes 

In this section we follow the previous analysis [Sj and 
provide a similar semi-analytic expression for the grav- 
itational wave strain emitted in the final stage of the 
spindown. An analytic result for the gravitational wave 
strain, averaged over polarizations and the orientation 
and position of a source at distance D on the sky, was 
given in ^ i38j 



h(t) 



3 oj^S22 



SOtt D 

with the mass current quadrupole moment 



5*22 = V2 — GManR^J 
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(36) 



(37) 



However, those authors did not determine the cutoff fre- 
quency analytically. Using our result eq. ( 26 1 we find 
for the final gravitational wave strain emitted in the last 
stage of the evolution 



3553 zT 



/ 



(38) 



in terms of the final frequency il/ eq. ( 26 ) . For a neutron 
star the final gravitational strain is 



.p.7.23.10-al(ME^) 



(39) 



The emitted gravitational wave strain during the evolu- 
tion is shown in fig.[5]for different values of the saturation 
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Figure 5: The gravitational wave strain of a 1.4 APR 
neutron star located somewhere in the Virgo cluster (distance 
20 Mpc). Shown are the numerical curves for two different r- 
mode saturation amplitudes ctsat = 1 (top) and Qsat ~ lO"** 
(bottom). The dots denote the analytic result eq. (381 for 



the minimum strain amplitude of the gravitational radiation 
emitted at the end of the r-mode evolution and the endpoints 
for other amplitudes lie on the line through these two points. 



amplitude. The numerical curves are compared to the an- 
alytic results for the minimum gravitational strain at the 



end of the evolution eq. (39) denoted by the points. As 



can be seen small amplitude r-modes have a lower max- 
imum amplitude but emit low amplitude gravitational 
waves over longer times. 

For an actual detection a systematic search strategy is 
required which takes into account the periodic structure 
of the signal. In ref. [S] matched filtering was consid- 
ered as a possible method. Since it requires prior knowl- 
edge of the detailed form and time evolution of the signal 
which might not be feasible, the estimates for the signal- 
to-noise ratio using this method are an optimistic esti- 
mate for more realistic methods. To compare the grav- 
itational strain to current detector bounds the Fourier- 
transformed gravitational wave strain h (i/) in the fre- 
quency domain is the relevant quantity. It can be ob- 
tained in a stationary phase approximation as 



dt 
dv 



(40) 



As seen from fig. [s] there are two qualitatively differ- 
ent stages of the gravitational wave emission. The first 
arises from the growth phase of the r-mode and results in 
narrow spike in h{i') at the stars initial frequency. This 



phase lasts only for a short time of order minutes and as 
shown in [8j there is not enough energy in this pulse to be 
detectable. Therefore, we neglect this initial part of the 
spectrum here. In the subsequent saturated phase where 
the star spins down it emits gravitational waves over a 
continuous frequency range with corresponding strain 



l9IGMR^ 



(41) 



The quantity that can directly be compared to the detec- 
to r sensit ivity, described by the rms strain noise hims = 
vSh{i^) in terms of the power spectral density^ Sh of 
the detector noise, is the characteristic amplitude of the 
signal defined by he = i' h. In particular the minimum 
characteristic amplitude hcj = h^Vf) reached at the 
lower frequency boundary of the evolution reads in the 
saturated regime 



T-c.f 



' 3IGMR^ 




+<y-2l3 



(42) 



^ ^^QCD 



in 

1/(4(39+4(7+7(7-60)) 



The square root in eq. (41) makes the dependence of 
he J on the microscopic parameters even weaker than for 
the final frequency eq. (29). However, the expression de- 



pends more strongly on macroscopic quantities like the 
mass and the moment of inerti a en coded in the parameter 
/. The rigorous bounds eq. ( 13 1 show that this depen- 



dence is also moderate since all these quantities vary only 
within roughly a factor of two. For a standard neutron 
star we find 



3/46 



5/368 




j^ j X 63/184 / ^ X 281/368 

1.4 Moj Vll.Skmj 



-29/184 



/Mpc\ _5/i84 



Note that this quantity is nearly independent of the satu- 
ration amplitude Ogat (and even increases with decreasing 
amplitude). The same holds for the microscopic parame- 
ters. For example the dependence on the neutrino emis- 
sivity involves the power 5/368 which provides a striking 
example of how insensitive such macroscopic observables 
can be to microscopically decisive yet inherently poorly 



® The "power spectral density" is just the square of the amplitude 
spectral density of the noise expressed in equivalent gravitational 
wave strain, as given e.g. by the LIGO collaboration [39|. 
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Figure 6: The characteristic gravitational wave amplitude of 
neutron stars, located somewhere in the Virgo cluster (dis- 
tance 20 Mpc) and based on the different APR stars given in 
tab.[llj compared to the sensitivity of the LIGO (dotted) and 
advanced LIGO (dashed) [9] detectors. The large dots denote 
the analytic expressions for the lower frequency limit for an 
amplitude Osat = 1 and the smaller ones for lower amplitudes 

down to Oaat = 10^''. 



known quantities, since even a drastic change by three 
orders of magnitude would only have an irrelevant effect 
of less than 10%. In order for the characteristic ampli- 
tude to give a meaningful estimate the observation time 
should be comparable to the time scale of the frequency 
change of the signal (so that all strength of a given fre- 
quency is actually detected) . Clearly for small saturation 
amplitudes where the spindown time increases drastically 
with decreasing ctgat i see eq. ( 35 1 and table III this is not 



guaranteed and therefore the fact that the characteristic 
amplitude is above the background alone is not sufficient 
for a detection. 

The result for the characteristic amplitude in the satu- 
rated regime is given for different stars in fig.[6j Whereas 
the spectra of the two standard neutron stars with modi- 
fied Urea reactions have nearly the same lower frequency 
cutoff, the heavy star with direct Urea interactions in 
the core has a higher cutoff > 100 Hz, but otherwise fea- 
tures a characteristic amplitude of similar size, as could 
have been expected from the insensitivity to the neutrino 
emissivity in eq. ( 42 1 . 



B. Theoretical signal detectability 

The signal to noise ratio for matched filtering can be 
written as |8j 



(tot.; 



ax 



91 
10L»2 



dv 



(43) 

in terms of the power spectral density St of the noise 
and ranging over a range from the final frequency of the r- 
mode evolution Vmin = i^f to the initial frequency v^ax = 
Vi at which the star is created. Here the label (tot.) has 
been added to indicate that the total frequency spectrum 
is included in this quantity. This point will be discussed 
in more detail below. In eq. (43 1 logarithmic frequency 



intervals are weighted by 1/5/,., since at low frequencies 
the frequency changes slowly which makes it easier to 
detect the periodic signal. Based on early estimates of 
the noise strain of a next generation LIGO detector and 
the assumption that r-modes are absent below 10^ K, 
resulting in a large Vmin ~ 120 Hz, the initial analysis 
[S] approximated the background noise by simple power 



law Sv, 



Consequently, they found that the integral 



is completely dominated by the low frequency end and 
the dominant quantity is the minimum frequency - which 
would surely be most welcome in view of the analytic 



expression eq. ( 29 1 . They find for the anticipated second 



generation LIGO detector ^ 



S_ 
N 



(tot.; 



Owen ct. al. 



1/2 



/ D 



lO^Scgsy V20Mpc 



V120Hzy 
(44) 

Yet, the anticipated design spectrum of the advanced 
LIGO background jH] improved compared to these early 
estimates. The noise has a minimum at v 300 Hz 
and increases again towards lower frequencies. There- 
fore, even for a large lower cutoff frequency assumed in 
[8J the above integral is now dominated by intermedi- 
ate frequencies which requires a more careful evaluation. 
After a numerical integration an approximate fit to the 
result over the relevant parameter ranges for I'min Eind 
i^max is given by^" 



S 
N 



(tot.) 



4.09 • lO^^vlfe- 



aLIGO 
X ^1+ 0.40 log 



D 



1000 Hz 



0.36 



(45) 

100 Hz - Vmin 

100 Hz 



i.e. we find a much weaker dependence on the minimum 
final frequency and a stronger dependence on the un- 



We fit the integrated form rather than the initial strain noise S^, 
since this provides a more precise result. Fitting the strain noise 
with a double power law Sh (i^) ~ au~^ + bv^ + c, similar to 
the ansatz used in |35| before their subsequent approximation, 
we can perform the integration analytically, yielding the same 
qualitative leading parameter dependence on Vm in (power law) 
and Vmax (logarithmic) we use for the fit in eq. (|45|. 
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known maximum initial frequency. By inserting the re- an average frequency v. The change l^v <^ v during this 



suit for the minimum frequency eq. (29 1 we see that the time interval can be obtained from eq. (19) 



uncertainty of this result due to the microphysics is mi- 
nor. The main uncertainty arises from the moment of 
inertia which using the above bounds has a maximum 
theoretical uncertainty of a factor of 4 but the actual 
range should again be smaller. 

The signal to noise ratios for different sources at differ- 
ent distances are shown in table |IV| Here we compare the 
values for sources given by the different considered stars, 
namely the IAMq and 2Mq APR stars for which only 
modified Urea reactions are possible and the 2.2 Mq star 
where direct Urea reactions are present within its core. 
These are assumed to be located at different distances, 
namely in the Virgo Cluster (20 Mpc), the Local group 
of galaxies (1 Mpc) and within the Milky way (30 kpc). 
Using the actual anticipated noise background of the ad- 
vanced LIGO detector and without the assumption that 
the evolution stops at 10^ K, our values are slightly larger 
those obtained in [8]. Heavier stars tend to give larger 
values due to the larger moment of inertia whereas di- 
rect Urea processes reduce the signal to noise ratio since 
the thermal steady-state curve is at lower temperatures 
and its intersection with the boundary of the instability 
region is correspondingly at a higher frequency that sets 



the lower gravitational wave frequency cutoff in eq. ( 45 I 



Yet, overall these values are rather similar for the differ- 
ent considered stars. These considerations suggest that 
theoretically a detection is possible for sources up to large 
distances in case the gravitational wave signal could be 
detected over its complete frequency range. 



C. Practical signal detectability 

Although the above signal to noise ratios might give 
some measure for the relation of the total signal that 
is observable from such sources and the detector back- 
ground, they do not directly reflect the detectability of 
the signal of a given single source over the actual ob- 
servation period, as noted already in [5]. They would 
provide such an actual measure for the detectability only 
for a short multi- wavelength signal that is completely ob- 
served, for example a binary inspiral. Yet, they cannot 
give a reasonable measure in the present case where the 
source is monochromatic and changes slowly with time 
while it is only observed for a short time interval, since 



the frequency ranges in eq. (43) that would only be de 



tectable hundreds of years before or after the present ob- 
servation can hardly affect the detectability within the 
limited observation interval. Therefore a proper signal 
to noise ratio that measures the detectability must take 
into account the finite duration of the observation and 



only include the relevant frequencies in eq. ( 43 1 



To find such a reliable signal to noise ratio for the 
detectability in gravitational wave detectors, assume that 
the frequency of the source changes only slightly during 
the short observation time interval Ai over which it has 



37r dt 



tg 



n) 



At 



(46) 



The current angular velocity Cl can be obtained in terms 
of the current age t of the star from the solution of the 
spindown equation eq. ( 19 ) which yields in the limit fli 3> 



Av 



V At 
6T 



(47) 



For Av <^ V the integral in eq. (43 1 gives an estimate for 



the actual observable signal to noise ratio 



S_ 
N 



(obs.) 



IqIGMR'^ Av 



SGI 



At 



(48) 



This quantity depends on the observation time interval 
At (see also ^40j) and depends on the internal composi- 
tion of the star only via the moment of inertia with the 
resulting uncertainty for the signal to noise ratio of a fac- 
tor of 4 noted above. Moreover, it depends only on the 
present (averaged) gravitational wave frequency v and 
not the complete spectrum. Strikingly it is again not 
directly dependent on the saturation amplitude which 
only enters indirectly via the dependence on the current 
frequency via Sh{v). Yet, we stress again that this ex- 
pression is only valid during the later evolution where 
the spindown becomes sufficiently slow so that Av <^ v. 
The right part of table |IV| shows the values for the ad- 
vanced LIGO detector for the different stars and a ratio 
At/t = 10~^ compared to the previous signal-to-noise es- 
timates. As can be seen in this case the observable signal 
to noise ratios are more than two orders of magnitude 
smaller than the total expressions eq. ( [45| . 

For the advanced LIGO detector the background is 
nearly constant over the relevant range from 80 to 800 
Hz with Sf/^ ~ 4 • 10~^^ and thereby allows us to further 
approximate to obtain the simple estimate 



S 

N 



(obs.) 



aLIGO 



lAt / 20 Mpc 



(49) 



I.e. for a one year observation period, comparable to past 
LIGO runs, a source in the Virgo cluster (20 Mpc) would 
be impossible to detect unless it is very young. A source 
in our local group of galaxies (1 Mpc) in contrast might 
be detectable if the results obtained assuming matched 
filtering do not strongly overestimate the detectability 
for more realistic search methods. A source in the Milky 
Way (30 kpc) in contrast should be detectable even if a 
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neutron star 


Tf[K] 


// [Hz] 


t^f [Hz] 


tsd [y] 


hf @ 1 Mpc 


hcj @ 1 Mpc 


iV 1 Virgo 


JvIl.r. 


c 1 1 tni 1 

iV Im. w. 


S 1 

iV 1 VirRO 


S 1 {.o^^) 


S 1 (o'^s) 
1 m. w. 


APR 1.4 M© 


1.26 ■ 10^ 


61.4 


81.8 


12.3 


7.2 ■ lO"^'^ 


1.6 ■ 10"^^ 


11.7 


234 


7794 


0.084 


1.7 


56 


APR 2.0 M© 


1.39 • 10^ 


60.3 


80.4 


4.9 


9.7- lO"^'^ 


1.8 • 10~^^ 


14.0 


281 


9356 


0.099 


2.0 


66 


APR 2.21 Mq 


3.15 ■ 10** 


84.0 


112 


1.7 


2.2 • lO"^'' 


2.1 ■ 10"^^ 


12.7 


254 


8460 


0.094 


1.9 


62 



Table IV: Results for the spindown and gravitational wave observables for different neutron stars considered in this work 
obtained from the semi-analytic expressions eqs. \27\ , ( |26[ ), ( |34[ |, ( |38[ ), ( |42[ ) and ( |43[ ). In contrast to the lower mass stars, the 
2.21 Mq neutron star is dense enough to allow for direct Urea reactions. All values are given for a large saturation amplitude 
otsat = 1. The signal to noise ratios are given for the advanced LIGO detector and are obtained with the gravitational wave 
frequency corresponding to the Kepler frequency as an upper cutoff of the spectrum respectively a fiducial observation interval 
Ai/t = 10"^ They are given for sources at three different distances ranging from the Virgo cluster (20 Mpc), and the local 
group of galaxies (1 Mpc) to a source within the Milky Way (30 kpc). 



more realistic search strategy has a significantly lower sig- 
nal to noise ratio. Since the original LIGO detector had 
a more than an order of magnitude larger background it 
is less surprising than for the previous estimate eq. (44) 
that no signal has been detected so far. As shown in 
the previous section the r-mode spindown scenario re- 
quires a sufficiently large amplitude leading to a suffi- 
ciently short spindown time so that after a long obser- 
vation time of several years, the sensitivity could still be 
sufficient to detect pulsars that are currently in the r- 
mode phase. Future third generation gravitational wave 
detectors like the currently planned Einstein telescope 
|41| are expected to exceed advanced LIGO in sensitivity 
by an order of magnitude and thereby should allow us 
to see many more potential sources. This should mark 
the advent of gravitational wave astronomy and could al- 
low us to learn about the composition of compact objects 
by gravitational wave observations. In this endeavor the 
simple semi-analytic expressions given here might prove 
useful. 



Finally let us estimate likely values for the above sig- 
nal to noise ratio for a few examples of known possible 
sources that could be young enough to emit gravitational 
waves due to r-modes. These consist of, in addition to 
the fastest pulsar J0537-6910 discussed above, the young 
neutron star Cassiopeia A from which no pulsation has 
been observed and several supernova remnants where no 
compact object has been detected so far. However, it 
has recently been suggested that the hard x-ray emission 
from SN 1957D results from a pulsar wind nebula and the 
compact remnant which produced it would then be the 
youngest indirectly observed pulsar [42]. The bounds for 
the moment of inertia again introduce an uncertainty of 
at most a factor of 4. As can be seen despite the fact that 
the observable signal-to-noise ratio eq. ( 49 1 is reduced by 



more than an order of magnitude compared to the to- 
tal versions eqs. (|44| and (45 1 the values are still large 



since all of these object are rather near. Even though the 
original LIGO had an order of magnitude weaker sensi- 
tivity the values are with the exception of SN 1957D and 
J0537-6910 large enough that it should have seen gravi- 
tational wave emission from these objects. Therefore it is 
very likely that these objects are not compact stars that 
spin fast enough to be in the unstable region. The proba- 



bly youngest pulsar, produced by SN 1957D, would be an 
ideal candidate, but it is most likely too distant to clearly 
detect gravitational waves due to r-mode emission with 
advanced LIGO in its standard configuration. However, 
the further improved sensitivity of advanced LIGO in a 
narrowband mode [S] could increase the signal-to-noise 
ratio above the detection level. Therefore, focussing on 
the frequency range in table [V] this source might present 
a promising candidate for such an in-depth search, pro- 
vided that the observation interval is long enough that 
the sensitivity gain is not outweighed by the limited data. 
In case the fastest directly observed young pulsar J0537- 
6910 is still in the final stage of its r-mode evolution there 
is a good chance that its signal was just below the thresh- 
old of the original LIGO detector and will be detected by 
advanced LIGO. 



VII. CONCLUSIONS 

We have derived general semi-analytic results for the 
key quantities of the spindown evolution of young pul- 
sars that reveal the dependence of these observables on 
the relevant physical ingredients. Vvfe find that the final 



spin frequency eq. ( 26 ) and the minimum gravitational 



wave strain emitted at the end of the evolution eq. ( 38 1 



are remarkably insensitive to the microphysical material 
properties. For other input quantities like the moment of 
inertia of the star or the angular momentum of the mode 
we gave rigorous bounds that allow us to estimate the un- 
certainty in these parameters. We show that the thermal 
evolution is always faster than the spindown evolution. 
Therefore, if the form of the instability region is convex 
a steady state is reached where viscous heating equals 
cooling due to neutrino emission. Non-convexity of the 
instability region due to stability windows arising from 
enhanced dissipation due to exotic phases of matter can 
prevent the star from reaching the thermal steady state 
and thereby drastically change the evolution which could 
lead to clear astrophysical signatures. 

The semi-analytic results for the final frequency of the 
evolution enabled a quantitative comparison to pulsar 
timing data and we find that all stars except one lie 
clearly below the frequency to which r-modes can spin 
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source 


age [y] 


D [kpc] 


u 


h 


he 


s (of's) 
N LIGO 


N aLIGO 


s ('°*) 

-AT aLIOO 


J0537-6910 


7 


52 


83 


1.4(0.1) X 10"^® 


6.5(6.8) X 10"^° 


2.1(0.2) 


28 (2.8) 


- 


SN 1987A 


25 


51.5 


157 (337) 


9.9(1.0) X 10"^^ 


0.8(1.2) X 10"^^ 


23(18) 


221 (236) 


2.0(3.2) X 10=* 


SN 1957D 


55 


4610 


137(296) 


0.7(2.3) X 10"^** 


0.9(1.3) X 10"^^ 


0.2(0.1) 


1.6(1.8) 


19 (33) 


Gl.9+0.3 


140 


7.7 


117(157) 


2.8(0.3) X lO-^'^ 


4.9(7.2) X 10"^^ 


57(62) 


598 (662) 


0.9(1.8) X 10'' 


Cassiopeia A 


« 300 


3.4 


104(223) 


4.3(0.4) X lO"^'' 


1.0(1.5) X 10"^* 


83 (106) 


906(1015) 


1.3(3.7) X 10* 


SN 1604 


400 


6.1 


99 (213) 


2.1 (0.2) X lO"^'^ 


5.7(8.3) X 10"^^ 


39 (53) 


433 (488) 


0.6(2.0) X 10'' 


SN 1572 


^ 440 


2.8 


97 (209) 


4.3(0.4) X 10"^^ 


1.2(1.8) X 10"^* 


80(111) 


894(1013) 


1.0(4.3) X 10" 



Table V: Expected gravitational wave signal for different young sources in case there is an associated fast rotating compact 
object that spins down via the r-mode mechanism. All values are given for the 1.4 Af© APR star model with Osat ~ 0.1 (0.01) 
- for Osat = 1 r-modes would have had already decayed at the current age of these systems. The observable signal-to-noise 
ratios are computed assuming an observation time At = 1 y. The frequencies z/ are obtained from eq. ( 33 1 in the limit that the 
unknown initial frequency was much larger than the current one. The theoretical total signal-to-noise ratios are here obtained 
by taking into account the total frequency interval [vf, u] that is observable from the stars current age to the end of the r-mode 
spindown. 



down a young pulsar. Only the fastest spinning young 
star PSR J0537-6910 is fast enough that it could lie 
within the instability region. Using additional observa- 
tion of spindown rates we find that this is unlikely since it 
lies just at the lower border of the uncertainty range. The 
spindown time in contrast depends strongly on the satu- 
ration amplitude and increases quickly from time scales 
of order years at asat ~ 1 to more than a million years 
below asat < 10~^. 

The r-mode scenario can therefore provide a quantita- 
tive explanation of the low spin frequencies of young pul- 
sars if the saturation amplitude is within certain limits. 
On the one hand it has to be sufficiently large asat ^ 0.01 
for the spindown to be fast enough to dominate other pos- 
sible spindown mechanisms and result in spindown times 
that are shorter than the ages of observed pulsars of more 
than several hundred years. On the other hand, even tak- 
ing into account the significant uncertainties, a very large 
saturation amplitude asat ~ 1 would have spun PSR 
J0537-6910 down to its current frequency in less than 
a hundred years. This would be in clear conflict with 
the much larger age estimate of its remnant of roughly 
5000 years [351 and the significant present spindown rate 
of the pulsar, which is most likely due to the electromag- 
netic radiation and should have further sped up the spin- 
down. Different saturation scenarios do not significantly 
change the spindown time and therefore should not qual- 
itatively affect our bounds on the saturation amplitude. 
In case r-modes reach sizable amplitudes asat ^0.1, PSR 
J0537-6910 could be younger than its characteristic age 
- obtained within a pure magnetic dipole model - sug- 
gests. If in contrast r-modes saturate already at ampli- 
tudes ttsat < 0.01 they would be subleading compared to 
the other spindown mechanisms, which must be sizable 
according to the significant spindown rates of observed 
young pulsars. 

Since we find that likely no observed pulsar is at 
present spinning down due to r-modes, one has to be 
careful when interpreting spindown limits on r-mode am- 
plitudes from young stars ^35, . Although these limits give 



important information on these stars, they are trivially 
fulfilled if the r-mode has already decayed. In this case 
the current rate involves only information on other spin- 
down sources, like magnetic dipole emission, and cannot 
give a bound on the r-mode saturation amplitude dur- 
ing the era when r-modes were actually present. The 
spindown limit [35] simply extrapolates the solid lines in 
fig. |4] outside the instability region. In reality for a large 
saturation amplitude the star would actually spin down 
quickly inside but then the spindown rate would drop to a 
lower value determined by the other sources once the star 
leaves the instability region. The slope of the curves for 
other spindown mechanisms are generally smaller, as il- 
lustrated in fig. |4] for the Crab and the Vela pulsar where 
a reliable braking index is available, and can therefore 
be perfectly consistent with observations even for a large 
saturation amplitude during the r-mode phase as long as 
the star is currently already sufficiently far away from the 
instability boundary. Yet, in case of J0537-6910 which is 
still close to the instability boundary, the spindown limit 
is more restrictive and coincides with the limits above 

m- 

We considered in this work only the dominant to = 2 r- 
mode. The inclusion of higher multipoles might slightly 
affect the spindown time by a factor < 1 but will not 
change the final frequency since the corresponding insta- 
bility regions of these modes are smaller |3] and there- 
fore the final segment of the steady-state curve, where 
the evolution spends by far the longest time, is identi- 
cal. Furthermore, we have here mainly studied the sim- 
plest case of a neutron star with minimal damping due to 
leptonic interactions and modified Urea processes in the 
core. The main reason for this is that it is the case of min- 
imal damping whereas other phases of matter or struc- 
tural inhomogeneities will increase the damping. There- 
fore, the final spindown frequencies we find are the lower 
limit to which general compact stars can be spun down 
by r-modes. In reality several aspects like pairing, mu- 
tual friction, the crust and other spindown mechanisms 
complicate the analysis. Although we did not study these 
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in detail here, our general semi-analytic expressions allow 
us to study any single source of damping, like the effect 
of an Ekman layer. 

Nuclear pairing could strongly change the temperature 
dependence of material properties over the considered 
temperature regime so that it would not have a simple 
power law form. Therefore it should alter thermal evolu- 
tion, but since at saturation the thermal evolution does - 
aside from the fact that its final frequency at the bound- 
ary of the instability region depends on temperature - not 
affect the spindown, the results presented here should not 
be strongly affected by pairing. The relevant part of the 
boundary is determined by shear viscosity which is dom- 
inated by electromagnetic leptonic processes which are 
only indirectly affected by the presence of hadronic mat- 
ter. As found in [23J this is a moderate effect and taking 
into account the insensitivity of our results to the shear 
viscosity this should only slightly increase the considered 
error range. 

As discussed before mutual friction in superfluid sys- 
tems is an additional source of damping that could con- 
siderably shrink the instability region. If this effect is 
very strong it could significantly raise the frequency to 
which r-modes can spin down a star. As can be seen from 
fig. [l] if the enhanced damping is only present at low tem- 
peratures, e.g. T < 10^ K as assumed in [5^, this would 
strongly affect the evolution at low saturation amplitudes 
but only moderately at sufficiently large amplitudes fa- 
vored above. 

Furthermore, it has been argued that the boundary 
layer at the interface between the core and the crust 
strongly enhances the viscous damping which would re- 
duce the r-mode instability region [43^ . This mechanism 
is based on the picture that there is a sharp boundary 
and the fluid velocity drops to zero over a distance of a 
few centimeters. However, in reality the crust is a com- 
plicated structure with possible geometrically structured 
pasta phases and an inner crust where the concentration 
of the neutron fluid decreases continuously. The tran- 
sition from a pure fluid to a rigid lattice then happens 
gradually over a large distance and it is questionable why 
the above picture should apply here. Phenomenologically 
these effects have partly been modeled via a slippage pa- 
rameter and the impact becomes rather small for likely 



parameter values [33]. We note, however, that our gen- 
eral analytic expressions should also apply when the crust 
provides the dominant source of dissipation. In this case 
further complications like the breaking or melting of the 
crust could be relevant which in turn should reduce the 
effect of viscous boundary layer damping. 

Finally, we computed the gravitational wave emission 
due to the r-mode spindown of young stars. We updated 
the expressions previously given in |S] by employing the 
current advanced LIGO design background and flnd that 
this further decreases the dependence on the flnal fre- 
quency of the spindown analysis and thereby on the mi- 
crophysics. Beyond this, we point out that a realistic 
signal-to-noise ratio should take into account that the 
signal is only observed for a limited time. The corre- 
sponding novel expression for the signal-to-noise ratio eq. 
( 48 1 is proportional to the root of the ratio of the obser- 



vation time and the age of the star and thereby becomes 
small for short observation times or old stars. For typi- 
cal values we find that the signal-to noise ratio is roughly 
two orders of magnitude smaller than the previous result. 
Although advanced LIGO would still be able to detect 
known nearby sources if they emit gravitational waves at 
present, this reduction limits the distance to which we 
can see such objects and thereby the number of observ- 
able sources. Nevertheless, both the steadily growing list 
of known young pulsars and the strongly increased sen- 
sitivity of next generation gravitational wave detectors 
like advanced LIGO could still make it possible to detect 
gravitational radiation due to r-modes, which would pro- 
vide a unique way to learn about the internal composition 
of compact stars. 
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